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I. INTRODUCTION 

Classical field theory in 1+1 dimensions has been used 
recently as a test ground for various approximation meth- 
ods for calculating the dynamical evolution of a nonequi- 
librium system ||l|, ^, ^. In classical field theory one 
can numerically simulate the full time evolution of the 
system via Monte Carlo methods Given the proba- 
bility distributions in the initial state of the system, one 
can randomly sample the initial distributions and numer- 
ically solve the equations of motion. The time evolution 
of the various expectation values is then obtained by av- 
eraging over the number of independent realizations. The 
dynamics becomes exact in the limit when the number 
of initial samples goes to infinity. 

It is important to be able to carry out calculations 
both in coordinate (lattice) and momentum space. In a 
lattice calculation, operators are numerically realized by 
using finite-difference representations of derivatives and 
integrals. Finite-difference methods, though leading to 
sparse matrices, are notoriously slow to converge. Thus 
the need of using higher order methods, like the nonuni- 
form grid Chebyshev polynomial methods |Q, ^ we used 
in the past few years, which belong to the class of spec- 
tral methods. Then the resulting matrices are less sparse, 
but what we apparently loose in storage requirements, we 
gain in speed. We do in fact keep the storage needs mod- 
erate also, as we can achieve very good accuracy with a 
moderate number of grid sites. Spectral methods have 
made possible for the first time to carry out complex dy- 
namical calculations at next to leading order in quantum 
mechanics H, ^, 0] and field theory ||] . 

In order to evaluate the quality of the various ap- 
proximations, one needs to find ways of comparing the 
truncated dynamics predicted by a given approximation 
scheme with the exact result. Hence the importance of 
classical field theory where Monte Carlo methods allow 
us to probe the true dynamics of a hamiltonian system. 



In practice, the exact calculation is always carried out in 
two steps: First, we consider a truncated space, which 
is characterized by a cutoff in momentum space. Once a 
solution is obtained for a given finite value of the cutoff, 
we increase the value of the cutoff and repeat the calcu- 
lation until the result becomes independent of the cutoff. 
This second step is what we usually refer to as taking the 
continuum limit. 

Lattice and momentum space (continuum) calculations 
do give the same result in the continuum limit. However, 
in a truncated space the results will be different. Since 
we are trying to compare the exact dynamics of a system 
with the result of a non-lattice based calculation, we can- 
not reliably use the lattice result. Hence the need for a 
continuum implementation of the Monte Carlo method. 

It is interesting to remark here that the study of the 
classical (j)'^ field theory problem represents a warmup for 
the next-to-leading order studies of the quantum linear 
sigma model and beyond. For the linear sigma model, 
the presence of the Landau pole will prevent us from 
taking the momentum cutoff to arbitrarily large values. 
Therefore, we are in fact less interested in the continuum 
limit of the problem, at least for now. 



II. THE LAGRANGIAN 

For a real field ^(x), the classical Lagrangian density 
in the (f>'^{x) theory is given by: 

^('^,9^'^) = i [a^0(x)][a^</>(x)] - ^ [cl^^x) - rl]' , 



£(0, 9^0) 
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where ~ — ArQ/2 > 0. The classical equation of mo- 
tion for (pix), in 1-f 1 dimensions, is: 
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d^~dl + fi' + ^cj,'{x,t) 



{x,t) =0. 



(2) 
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In order to make contact with the work of Aarts et 
al 0, in this paper we address only the symmetric case, 
(^(0)) = (7r(0)) = 0. 



III. INITIAL VALUES 



Assuming that initiaUy the system is in thermal equi- 
librium, the initial values 0(x,O) and 7r(a;,0) = 0(a;,O) 
are taken from a canonical ensemble governed by a clas- 
sical density distribution tt] defined as 



p[(/.,7r] = Z-i(/3o)e-^«^[^--l, 



(3) 



with /3o = I/Tqi and 7r(x,t) = (j){x,t). Correspondingly, 
the ensemble average of a quantity tt] is defined by 



{A[(f>,Tr])^TT{pA[^,7T]} 



(4) 



Following Aarts et al [Q, we choose initial values by 
randomly sampling the density distribution correspond- 
ing to the free particle Hamiltonian (A — 0). For each 
set of initial conditions, we time-evolve (j){x,t) using the 
equation of motion (^). The average value 4'1i{t) is cal- 
culated as 



Ma 



Elim 
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i—1 



(/)2(x,i)dx, (5) 



where (i) denotes the i^^ Monte-Carlo run and Mc is the 
total number of Monte-Carlo runs. 



IV. MOMENTUM SPACE APPROACH 

In the momentum space approach one introduces the 
Fourier transform of the field 4>{q,t) via 

One then obtains the classical equation of motion in mo- 
mentum space 



The equation of motion is solved by using an Euler 
method, where the time-differential operator is replaced 
by a second order difference formula, i.e. 

d^4>{q, t) ^ t + T)~2 4>iq, t) + 4>{q, i - t)1 / r2 , 



where r is the time step. A momentum space cut-off is 
used and the system is also in a box of length L. Letting 
L — Na, and choosing our momentum to be 
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we then have that 



k-< 2'---'2 



1 



A = 



ttN 



(6) 



(7) 



where A is the momentum cutoff and a is the lattice 
spacing. Then, we must take the continuum limit as 



hm 



J -A 27r^' 



In this approach, the periodic boundary assumption is 
implied whenever we perform numerical Fourier trans- 
forms and convolutions, the continuum Fourier being re- 
placed by the discrete one. 

We will choose the initial values for (l){q, 0) and 7r(g, 0) 
to be solutions of the unperturbed Hamiltonian 



Ho 



dx I n^{x) + [d^Mx]] +M^0o(a;) 



(8) 



where ttqIx, t) — (j>Q{x,t). The corresponding equation of 
motion in momentum space is 



with the dispersion relation 

2 2,2 

General solutions of Eq. (^) are of the form: 
and 

TTo{q,t) ^ 4>Q{q,t) ^ age 



(9) 
(10) 
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We can now calculate the classical density matrix as 



^o(/3o) n' 



(11) 



■^o(/3o)=J| j j d<j)o{x) d7ro(a;)e" 



-/3o-Ho[0o,7ro] 



where we have put aq^ — Xq^ +Wqk ■ Hence, the symmetric 
case scenario, (0(0)) = ('''(O)) — 0, implies that Xq^ and 
j/qj. are uniform random deviates between 0.0 and 1.0 . 

Using Eq. with the density distribution given by 
(|rT|), we obtain 

{aq.aq^,) = (a,, a?,,) = 0- 

(K.o, ) = {aq^a ) = nq^{f3o) 6kM' / L , 
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where ng(/3o) = l/(/3o^q)- Note that nq{f3o) is the high 
temperature hmit of the classical Bose-Einstein occupa- 
tion number distribution. Finally, we obtain 



where 



1 1 



y - 



1 



/(^,A) 



1 1 _ 1 1 



with 



/(m,A) 



1 



arctan(A//i) . 



(12) 



Thus, for a given cut-off in momentum space, we have 

1 



(r(o)) = 



arctan(A/^) 



(13) 



/i/3o7r 

In the limit when A oo, we obtain {(fP'iO)) = l/(2/i/3o)- 

V. PERIODIC LATTICE APPROACH 

For the lattice formalism, we follow closely the ap- 
proach presented in Reference ||l| . We discretize the con- 
tinuum equations by using a lattice in coordinate space 
with spacing a and periodic boundary conditions. The 
differential operators are replaced by second order differ- 
ence formulas: 

dl(f>{x,t) ^ (l){x + a,t) - 2(f>{x,t) + (j){x ~ a,t) /a^, 

dt(t>{x,t) ^ <p{x,t + T)-2(l){x,t)+(j){x,t~T) /t^. 

The dispersion relation is modified due to the Laplacian 
on the lattice and becomes 



-2 I 2 

q +H , 



— (1 - cosag) , 



(14) 



with a = tt/A. The momentum q takes the same finite 
number of discrete values, see Eq. The relationship 
between the cutoffs in the two approaches is given by 



lim 



lim J- y 



lim 



.A 2^^ ^ 



The initial conditions are generated by sampling the 
initial probability distribution of the unperturbed sys- 
tem just as in the continuum case, with the formal dif- 
ference that we replace the values of the momenta by 
the shifted values g^. Then, the initial expectation value 
becomes 



2^1 + (/ia/2)2 
Thus, for a given lattice spacing a = tt/A, we obtain 



(15) 



M/3o 2yrT(Ma72? 



(16) 



In the limit when a ^ (A ^ oo), we recover the con- 
tinuum limit (0^(0)) = l/{2fi/3o). However, for a given 




FIG. 1: /(/i = 1, A) and I{p = 1, A) as a function of A. 

momentum space cutoff A, the values of {(j)'^{0)) given by 
Eqs. (|l3|) and (16) are not the same (see Fig. ^. This 
is an artifact of solving the equations of motions on the 
lattice and requiring 4>{x, t) to satisfy periodic boundary 
conditions. As a consequence, one cannot directly com- 
pare continuum and lattice calculations. 



VI. RESULTS 

We choose to illustrate the approaches presented 
above, for a set of parameters which allows us to com- 
pare with results available in the literature ||l| . We have 
A = 1/3, n = 1, A = in, To = 5.03891094. Then, 
the initial condition corresponding to the momentum 
space approach is (0^(0)) = 2.39208677, as obtained 
from Eq. (p3), while the numerically calculated initial 
average is 0^i(O) = 2.39187089, for a 0.009% error. In 
turn the lattice calculation produces the continuum limit 
((/)2(0)) = 2.5, provided that one choose the lattice spac- 
ing a = 0.25. The lattice spacing is subsequently left 
unchanged, even though one may vary the number of 
lattice sites, and implicitly the lattice size (see Eq. (0)). 

The numerical methods used for the numerical imple- 
mentation of the two methods are well under control. In 
Fig. |2| we show the Monte-Carlo value of (t), as defined 
in Eq. (^), obtained using the periodic lattice approach. 
The error lines represent coordinate average deviations of 
the runs as a function of t. Similar results are obtained 
using the momentum space approach. 
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FIG. 2: Monte-Carlo calculation of 4>ci{t) versus t. 




FIG. 3: Comparison of the lattice and the continuum Monte 
Carlo results, respectively. 



MC - lattice, shifted 
MC - continuum ■ 




VII. CONCLUSIONS 

In this paper we have discussed two approaches of ob- 
taining the dynamical evolution of a classical system, one 
based on a lattice formulation in coordinate space, the 
other in momentum space. Both methods require the 
assumption of periodic boundary conditions, but the dif- 
ferent levels at which this assumption is made, allows 
the momentum space approach to avoid certain artifacts 
of the lattice based method. In particular the intrin- 
sic mismatch in initial conditions at finite cut-off values, 
results in different values of the "thermalized" field, at 
large times. The discrepancy is worse for smaller values 
of the cut-off, but the two approaches converge to the 
same result in the continuum limit. The mismatch in 
initial conditions is due to the fact that by using a finite 
difference approximation for the spatial derivative oper- 
ator together with the assumption of periodic boundary 
conditions on the lattice, we have in fact introduced a 
approximation of the dispersion relation (see Eq. ( p^ - 
lattice, and Eq. ( p^ - continuum), which is now viewed 
as an expansion in the lattice spacing a. In order to im- 
prove the quality of the spatial derivative approximation 
in the lattice case, one would normally have to take the 
limit when the lattice spacing a goes to zero. We are 
however prevented from doing that, since the choice of 
the momentum cut-off A also determines the choice of 
the lattice spacing a — tt/A. Consequently we cannot 
improve the agreement of the lattice dispersion relation 
with the continuum for a given momentum space cut-off. 
The momentum space (continuum) approach does not 
exhibit this limitation. 

One may think of modifying the unperturbed Hamilto- 
nian in order to effectively obtain a higher-order approx- 
imation of the dispersion relation while still having the 
same type of equations of motion, similar to the improved 
action framework in lattice QCD. This would result in 
new values g^, and would require cancelling the various 
orders of a in a rigorous fashion. (The standard lattice 
calculation introduces values of which differ from 
already at order a^.) However, this is beyond the scope 
of the present work, and since we are in fact able to ob- 
tain an exact solution for the continuum problem, we 
merely state here the differences between the lattice and 
the continuum approach. 



FIG. 4: Comparison of the shifted lattice and the continuum 
Monte Carlo results, respectively. 
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